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Abstract 

We investigate the finite temperature behavior of the meson sector of an ef- 
fective Lagrangian which describes nuclear matter. A method is developed for 
evaluating the logarithmic terms in the effective potential which involves expan- 
sion and resummation; the result is written in terms of the exponential integral. 
In the absence of explicit chiral symmetry breaking, a phase transition restores the 
symmetry at a temperature of 190 MeV; when the pion has a mass the transition 
is smooth. At a much higher temperature a first order phase transition restores 
scale symmetry. 
PACS: ll.lO.Wx, 12.39.Fe 
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1 Introduction 



Spontaneous chiral symmetry breaking has long been studied in the hnear 
sigma model of Gell-Mann and Levy . At finite temperature the restoration 
of chiral symmetry in both the linear and non-linear versions of the model 
has been discussed in, for example, Refs. 0, |^ (see also references therein). 
If the potential of that model, ^A(cr^ + tt^ — /^)^, is used in calculations of 
nuclear matter it leads to compression moduli which are much larger than the 
observed value; furthermore the binding energies predicted for finite nuclei 
are much too small @, In previous work |^, hereinafter referred to 
as I and II respectively, we have therefore replaced this potential with a 
form which incorporates broken scale symmetry in addition to spontaneously 
broken chiral symmetry, as suggested by quantum chromodynamics (QCD). 
In particular the potential contains logarithmic terms involving the glueball 
field and the a and tt fields. At temperature T = this led to a good 
description of nuclear matter and finite nuclei at the mean field level. In II 
the Lagrangian was extended to include explicit chiral symmetry breaking, 
an additional chiral-invariant term, and the isotriplet vector mesons, and 
satisfactory agreement with low energy vrA^ scattering data was obtained. 

The purpose of the present paper is to examine the finite temperature, 
T > 0, properties of our Lagrangian. Previous studies 0, |] of models of 
this general type at T > have simply included temperature effects for the 
nucleons. Clearly thermal effects for the mesons are also needed, particularly 
those due to the pion which will be dominant at low temperatures. However, 
while the use of logarithmic potentials is straightforward at T = 0, it is far 



2 



from obvious how to proceed at T > 0, even at the mean field level. Since 
the analysis is quite complicated, we will focus here on the mesonic part of 
our Lagrangian which contains the 0, a and tt fields. The Lagrangian and 
our thermal analysis is discussed in Section 2. We give our numerical results 
in Section 3 and Section 4 contains our conclusions. 



2 Theory 

2.1 Equations of Motion 

As mentioned, we simplify the Lagrangian by excluding nucleons, as well as 
the oj meson which couples to them. We also take the simplest form for 
the mesonic contributions from 1, augmented by the explicit chiral symmetry 
breaking discussed in II so as to endow the pion with a mass. Then our 
effective Lagrangian involves the glueball field and the chiral partner fields 
a and tt and takes the form 



V 



V 00 ^ ^2 



V 



(7 + TT — 







1^-2 
ctq 



2(2 

(1) 



Here C = ^ iii the vacuum (j) = cpo, a = ao and tt — 0, regardless 
of whether or not the explicit symmetry breaking term e[ is present (an 
additional term, €3, was unfavored in II and is omitted here). Thus we have 
spontaneous, as well as explicit, chiral symmetry breaking. The quantities B 
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and 6 are parameters. For the latter, guided by the QCD beta function, we 
take 6 = 4/33 as in I and II. The logarithmic terms here contribute to the 
trace anomaly: in addition to the standard contribution from the glueball 
field 1^, there is also a contribution from the a field. Specifically the trace 
of the "improved" energy-momentum tensor is 



4e, 



(2) 



where $j runs over the scalar fields {0, a, tt} and the vacuum energy, evac = 
-lB<Ptil - 6) - e[. 

We take the vacuum glueball mass to be approximately 1.6 GeV in view 
of QCD sum rule estimates fill] of 1.5 GeV and recent lattice estimates |]12 



of 1.7 GeV. Since the mass is large in comparison to the temperatures of 
interest we shall neglect thermal effects for the glueball. We define the ratio 
of the mean field to the vacuum value to be x = 4>/4'o- Then Lagrange's 
equations for the glueball and a fields in infinite matter are: 



= 45ox'lnx-5o5x 



2x'ln 



V2 + 7r2^ 







'2a 













cr 



(T2 + 7r2 



(3) 



where we have defined Bq = 

We wish to take into account thermal effects for the a and tt fields. To 
that end we break a into a mean field part a and a fluctuation Act with mean 
value (Act) = 0. The mean value of the pion field (tt) is, of course, zero. We 
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write 



— (a^ + 2aAa + Aa^ + tt^) 



'0 



z/^ + 2z/Az/ + V , (4) 



where, as in I, u = a/ao and here Az/ = Aa/ao and i/?^ = (Acr^ + 7r^)/crQ. 
The treatment of the formal equations (^) at finite temperature is far from 
obvious. Simply expanding the fluctuations out to lowest order, as in Ref. 
[|T^ , will not properly treat cr^ + tt^ when it occurs in the denominator or in 
the logarithm in Eqs. (|^). (It is known that approximating the logarithm at 
zero temperature gives poor results.) We shall proceed in two steps. For the 
first step it is useful to bear in mind that at low temperatures z/ ~ 1 and the 
thermal average (1/?^) is small, while at high temperatures u is small but the 
thermal fluctuations are large. This means that the cross term IvAv of Eq. 

is small in both limits. Thus we expand out this term; odd powers can 
be dropped since the thermal average gives zero. This yields 



^ B.,x [-2X^ ln(.^ + + j^^, + j^^, + + 

+5ox'(41nx -6)- e[xi2u - - (V') - x') , 
1 2Az/2 Au^Au^ Su^Au^ \ 



= 5o5xV 



+Bo6x'iy-e[x\l-iy) . (5) 



Here the angle brackets indicate that we have taken a thermal average. For 
most purposes it is sufficient to truncate at 0{h'^ + i/?^)"^, however in the 
absence of explicit symmetry breaking, e[ = 0, there is a small region near the 
chiral phase transition where solutions cannot be obtained. This difficulty is 
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alleviated by going one power higher. The expansion parameter is of order 
4z/^(Az/^)/(z/^ + {'4>'^)Y which a posteriori we find to be < 0.05 - this is 
satisfactorily small. For the thermal average of the square of a fiuctuating 
field, for example for a component Tr^ of the pion field, one has the standard 
result 



Here /3 = 1/T is the inverse temperature and cj^ = fc^ + m*^, with m*^ 
the effective pion mass which will be discussed in the next subsection. In 
calculating thermodynamic integrals, such as this, we find it convenient to 



make use of the numerical approximation scheme of Ref. ||14 |. 

In order to evaluate Eq. (|^) we take the second step, which is to make a 
formal expansion in i/?^. It is sufficient to consider the logarithm in Eq. (|^) 
which gives 

.l)n+l /^2N"y 



(ln(z/^ + V')) = lnz/2 + (^ 



n=l ^ V 



oo 



In + + -!)!?/-". (7) 



n=l 

Here we have defined 



In obtaining Eq. we have used the counting factors described in the 
Appendix. These require that (Aa^) = (7r2)/3. This will be exact in the 
high temperature regime when chiral symmetry has been restored. At lower 
temperatures where we use the appropriate values for the thermal expectation 
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values in Eq. @, this will be approximate. However we do obtain the correct 
low temperature limit. We also remark that if the fluctuations in the a field 
are arbitarily ignored, the qualitative behavior is similar to the present case. 
Thus we believe the approximation is reasonable, although we would wish to 
have a more quantitative assessment of the errors involved. It would seem 
to be essential to take into account in some reasonable fashion the vertices 
in Eq. (J^) with large numbers of fields attached. 

While the series in Eq. (|^) is divergent, we regard it as a formal expansion 
which must be resummed before it can be evaluated. The counting factors 
for an odd number of field components lead to expressions involving the error 
function. For an even number of field components, four in the present case. 



the resummation requires the exponential integrals |jT5| 



oo , 



E^{y)=jdt-^. (9) 



1 



Matching the series of Eq. (|^ to the asymptotic expansion of the exponential 
integrals for large we obtain 

= In z/2 + (1 - y)eyEi{y) + 1 . (10) 

In the limit of low temperature [y oo) or high temperature {y — > 0) 
we obtain 



(ln(z/2 + ^2))^lnz.2 + Ml {y^oo) 

- ,Z 



2 , „/.2^\ , i„ ,,2 _^ _ 

MO.TGS^V'')) (1/--0). (11) 
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The low temperature expression simply corresponds to expanding the loga- 
rithm to lowest order, as it should. In the high temperature expression the 
numerical factor is 0.7631 = e^~^/2, where 7 is Euler's constant. 

The remaining quantities needed in Eq. (^) can be evaluated in analogous 
fashion or, more conveniently, by differentiating Eq. (|T0|) with respect to v"^. 
In order to consistently evaluate the counting we replace Ai/^ in Eq. (^) by 
ii/j^. The terms involving a fourth power of the fluctuations are evaluated 



in analogous fashion using the approximation 
/ \ / Az/V? \ / 



2\2 



Then Lagrange's equations can be written in final form: 



(12) 



= Bo6x' 



(-2 + 2y-y^ + 2y^ + + iy^) - 3 + y - fi/^ 



Bo5x^y 



{3y^ + Qy^ + y')eyE,{y) - 3 + y - 5y^ - y' 



(13) 



(14) 



Exact chiral symmetry restoration, z/ = 0, will only be obtained from Eq. 
(|14D when there is no explicit symmetry breaking, e[ = 0. We also note that 
Eqs. (|13D and ([I^ permit a solution x = corrsponding to scale restoration, 
in which case u is undefined. 
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2.2 Masses 

For each field we define the effective mass at finite temperature as the thermal 
average of the second derivative of the potential. This means that we only 
consider contributions arising from a single interaction vertex. Since the 
mixing between the glueball and the a meson is small, we neglect it here for 
simplicity. Specifically 

= \M^/ = + + \-^^ + (a2 + 7r2)2 / ' 

2 *2 /d'V\ , I BoSalx^ 2Bo6alx^\ 

={9^J = (BoS + e,)x +\-^^+ (,2 + ^2)2 ; ' 

<f>K = (^) = 45ox'(3 Inx + 1) + 3(e; - Bo6)x' + {Bo6 + e[y 



.2e[u + ( -6BoSx' In ( ) + W + e[)il^' ).(15) 



As we have remarked, we do not consider thermal fluctuations in the glueball 
field and so its mass does not enter the equations. However it will be useful 
to display the mass in Sec. 3. The a and vr masses are needed in evaluating 
("0^) as indicated in Eq. (|^). As a result the equations of motion and the 
expressions for the masses must be evaluated self-consistently. Treating Eq. 
(|T5|) as discussed in the previous subsection and using the equations of motion 
(|T3|) and ( [T^ ) for the case where u ^ and x 7^ we obtain 



BoSxy 

3z/ 



V 

4„,2 

2 



'l%y + 15y2 + 2y')e^E^{y) + 7 + ISy + 2y' 



a' v2 



^o^^tt' = . (17) 
= 4(5ox' + e» - 2{B^b + e;)(z/2 + (V^^^) . (18) 



It is straightforward to verify by taking the y ^ oo hmit that the zero 
temperature resuhs of II are obtained. In the case when z/, ?/ — ^ the a and 
vr masses become equal: 



a^mf almf = ^ ^ B,5x' " ^) ' (19) 



where for the last expression we have used Eq. ([T^ . In the case where there 
is no explicit symmetry breaking, e'^ = 0, at sufficiently high temperature 
1^ = and the masses are precisely given by the last expression in Eq. (|T9|). 

2.3 Thermodynamics 

The grand potential per unit volume can easily be written down: 

- = /V) + ^ ( dkk^ hn(l - e-^'^^) + 31n(l - e"^'^-)" 
V 2tx^ J L J 

-im:'(Aa^)-imf(7r^), (20) 

where uj^ = k'^ + m*^ and uj"^ = k"^ + m*^. The subtraction of the last two 
terms in Eq. (pO]) is necessary to avoid double counting |jT6|. The evaluation 



of the thermal average of the potential follows the discussion in Subsec. 2.1 
and we simply quote the result. 
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+l-\y + l,y^ + ly' + ^,y' + Inz/^J + \[B,{1 -5) + e\] . (21) 

We have added a constant term here so that (V) is zero in the vacuum. The 
pressure P is of course —VL/V. 

Now if one takes the partial derivative of Vl/V with respect to x or the 
equations of motion (p!3D and (|14D ought to be obtained. This is true if one 
ignores the dependence of the masses on these variables. If this is taken into 
account derivatives of the explicit m? terms in ( PO] ) cancel with derivatives 
of the Bose partition functions. Derivatives of (i/?^) do not cancel precisely. 
They would do so if all the terms arising from derivatives of the original 
logarithm of Eq. (|T]) were retained in the subsequent equations. However 
this gives additional contributions to the pion mass so that in the absence 
of explicit symmetry breaking the mass is no longer zero at low tempera- 
tures in violation of Goldstone's theorem. As it is necessary to approximate, 
one cannot have it both ways. We prefer to truncate the equations of mo- 
tion and the mass expressions at a given order and approximate the grand 
potential- we have verified that the additional terms needed to produce the 
exact equations of motion are small in comparison to the terms retained. 

In this spirit we ignore the derivatives of ("0^) in deriving the energy 
density which takes the simple form 

f = (V)-|mf(Aa^)-|mf(7r^) 

(22) 
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3 Results 



The explicit symmetry breaking parameter, €[, is chosen to yield the vacuum 
pion mass {e[ = (J^rn^) for one set of calculations and to be zero for another 
set. In order to fit nuclear matter saturation in these two cases we find Bq 
to be (334.9 MeV)^ and (342.6 MeV)"^, respectively. Based on the results of 
II we take (Tq = 110 MeV and 0o = 140.9 MeV. We stress that in evaluating 
(t/>^) the values of the sigma and pion masses from Eqs. (|TB|) and (|r^) are 
used. 

In Figs. 1 and 2 we display the sigma and glueball mean fields as a 
function of temperature. Here, and in the subsequent figures, the dashed 
line refers to e'^ > and the solid curve gives the = results. In Fig. 1 
the solid line shows a chiral phase transition and the dotted curve indicates 
a region of instability where the pressure is not maximized. Similarly the 
dotted and dash-dotted parts of the curves in Fig. 2 indicate regions of 
instability for the scale restoration phase transition. The pion and sigma 
masses are displayed in Figs. 3 and 4, respectively. 

Consider first the t'^ = results (solid curves) where u and m* become 
zero at the chiral phase transition temperature T^. Below the u = 
solution results in unphysical (imaginary) masses. The value of can be 
obtained by observing where the masses in Eq. (^) become zero, giving 
Tc = V^xcq- The effect of the glueball field is rather small here since x is only 
slightly less than unity at these temperatures (Fig. 2). With our parameters 
Tc = 187 MeV. Fig. 1 indicates a weakly first order phase transition, however 
we do not believe that our approximation scheme is sufficiently accurate to 
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determine the order of the transition. As one would expect, the behavior in 
the transition region is sensitive to the prescription used to define the masses. 
For example, we have remarked that if the equations are truncated at a lower 
order no solution is obtained for a small region in this vicinity. For T > 
the field u = and the pion and sigma masses are equal. 

There are numerous estimates of Tc in the literature. The standard chiral 
model yields Tc = v^/tt which with the pion decay constant = 93 MeV 
gives Tc = 132 MeV. Apart from the numerical factor, our value is larger 
because ctq has to be greater than in order to fit nuclei; see the discussion 
in I and II. We remark that, to a good approximation, the effect of modifying 
(Jo is simply to scale the temperatures on the abcissae of the figures. As with 
the standard model our a mass is zero at the critical temperature. Gerber 



and Leutwyler |T7] have studied two-fiavor quantum chromodynamics using 
an effective chiral Lagrangian to three loop order and estimate a value of 190 
MeV for with massless quarks. Finally we mention that the two-flavor 
lattice QCD results of Brown et al. [|1^] show a second order phase transition 
for massless quarks, but this is washed out by any flnite mass. The latter 
seems to be in line with our dashed curves with e'^ > where there is no phase 
transition and u smoothly decreases to a small value. Symmetry restoration 
takes place at somewhat higher temperatures when the pion mass is flnite, 
as would be expected. Figures 3 and 4 show that the pion mass smoothly 
increases from the vacuum value and the sigma mass initially drops, but then 
starts to increase at T ~ 350 MeV to become degenerate with the pion mass. 

There is a further interesting feature which occurs at high temperature, 
although we caution that the model may well not be reliable in this region. 
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Also, physically, deconfinement may have taken place since lattice calcu- 
lations suggest that chiral restoration and deconfinement occur at similar 
temperatures [1^ . Indeed with a crude lowest order treatment of gluons and 
massless quarks, taking the bag constant to be \e^s,c\, we find that the decon- 
fined phase is preferred for T > 170 MeV. In spite of these caveats we show 
(Fig. 2) that at a rather high temperature the solution with x ~ 1 becomes 
unstable and x drops to zero. We interpret this as a first order phase tran- 
sition which restores scale symmetry. With the present approximation the 
transition temperature is ~ 550 MeV. The physics behind this is indicated 
in Fig. 5 by a qualitative plot of Q/V as a function of x with e[ > (note 
that each curve is normalized to zero at x = 0)- At low temperatures, while 
there is a x = solution, the minimum corresponds to x ~ 1- At T = 530 
MeV the depth of these two minima become equal and beyond this the stable 
solution is X = 0. The unstable solution is indicated in Fig. 2 by the dotted 
(e'^ = 0) or dash-dotted (e'^ > 0) curve out to the point where the minimum 
and maximum of the potential, marked by dots in Fig. 5, coalesce. The x = 
regime should not be taken too seriously since all that remains of Q are the 
thermal partition functions for massless cr and tt mesons. Nevertheless we 
believe that a first order phase transition of this type is a general feature 
of the model. The critical temperature is rather high here because of the 



coupling of the glueball to the cr and vr mesons; Agasyan ||2^ has discussed 
the pure glue sector where the estimated critical temperature is lower. 

In Fig. 6 we plot the glueball mass as a function of temperature. (Note 
that when x = we have set u = since it is undefined and is small prior 
to scale restoration.) The glueball mass begins to drop significantly with 
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decreasing x fo^' T > 400 MeV, but remains much larger than m* and m*. 
Thermal fluctuations are not expected to be important, except, possibly, for 
the highest temperatures discussed here. 

Finally in Fig. 7 we show the pressure versus the energy density with the 
points labelled by the corresponding temperatures. At low temperatures our 
pressure is slightly negative since our approximation to the grand potential 
results in inaccuracies at the 1% level in the delicate cancellation of large 
terms. The chiral phase transition is just visible for the solid curve in Fig. 
7, and beyond this temperature interactions become less important and the 
massless gas result E/V = 3P is rapidly approached. 

4 Conclusions 

We have discussed the finite temperature behavior of the meson sector of 
an effective Lagrangian with which we have successfully described nuclear 
matter and finite nuclei. This is non-trivial because of the ln((j^ + 7r^) term in 
the effective potential. Our method of handling this term involved expansion 
and resummation of an infinite series with the final result cast in terms of 
the exponential integral. 

Our results showed that at sufficiently high temperature the mean value 
of the a field became small, signalling chiral restoration. Hitherto it has 
not been appreciated that chiral symmetry is regained in this type of model. 
In the absence of explicit chiral symmetry breaking a phase transition was 
obtained at Tc — 187 MeV; we do not believe that our calculations are 
sufficiently accurate to determine the order of this transition. In the physical 
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case where explicit chiral symmetry breaking was present and the pion had 
a vacuum mass a smooth restoration of chiral symmetry was found, the 
onset of which occurred at a somewhat higher temperature. Well beyond 
this a first order phase transition took place in which the mean glueball 
field dropped to zero, implying restoration of scale invariance. This is an 
interesting physical feature, but deconfinement is expected to have taken 
place before this temperature is reached and the application of our model at 
such a high temperature should be taken with a grain of salt. One would like 
to know how these results, particularly chiral restoration, are affected when 
nucleon degrees of freedom are present in addition to the mesons. This will 
be the subject of future work. 

We thank J. Schaffner-Bielich for useful comments. We acknowledge par- 
tial support from the Department of Energy under grants No. DE-FG02- 
87ER40328 and DE-FG02-94ER40823. G.W.C. thanks the University of 
Minnesota for a Doctoral Dissertation Fellowship. A grant for computing 
time from the Minnesota Supercomputer Institute is gratefully acknowledged. 



Appendix. Counting for a General Vertex 

We need to evaluate the thermal average of (V'^)", where 

N 

^'=j:^l (23) 

i=l 

which involves the sum over the squares of N fluctuating fields. This corre- 
sponds to evaluating diagrams of the type shown in Fig. 8 where, in general, 
the flower would have n petals. We assume that the thermal average {ijjf) 
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is independent of the label i, which amounts to assuming the masses of the 
particles involved are the same. Because of this assumption the result of 
taking the thermal average of each possible pair of fields at a general vertex 
can be written 

ii^'r) = Cr^i^Y , (24) 

where c„ is a number which gives the counting. Consider ('0^)"+^. The two 
additional fields can be averaged together. Alternatively we can break one of 
the n original pairs and combine the additional fields with them in 2 ways; 
this requires that the additional fields have the appropriate label i leading 
to a factor of 1/A^.Thus 



((V'T^^) = (i + f ) {^'){^') 



2^n 0^ 



Cn+1 = ■ (25) 



N 

Since Ci = 1, we obtain 
_ (iV + 2n-2)!! 

For the case at hand, AT = 4, = (n + l)!/2'*. 



References 

[1] M. Gell-Mann and M. Levy, Nuovo Cimento 16 (1960) 705. 
[2] C. Contreras and M. Loewe, Int. J. Mod. Phys. A5 (1990) 2297. 
[3] A. Bochkarev and J. Kapusta, Phys. Rev. D54 (1996) 4066. 

17 



[6: 
[7: 



[9 
[10 

[11 

[12 

[13 
[14 



E.K. Heide, S. Rudaz and P.J. Ellis, Nucl. Phys. A571 (1994) 713. 

R.J. Furnstahl and B.D. Serot, Phys. Rev. C47 (1993) 2338; Phys. 
Lett. B316 (1993) 12. 

G. Carter, P.J. Ellis and S. Rudaz, Nucl. Phys. A603 (1996) 367. 

G. Kalbermann, J.M. Eisenberg and B. Svetitsky, Nucl. Phys. A600 
(1996) 436. 

P.Papazoglou, J. Schaffner, S. Schramm, D. Zschiesche, H. Stocker and 
W. Greiner, |nucl-th/9609"035 . 



J. Schechter, Phys. Rev. D21 (1980) 3393; A.A. Migdal and M.A. 
Shifman, Phys. Lett. B114 (1982) 445. 

H. Gomm and J. Schechter, Phys. Lett. B158 (1985) 449. 

M.A. Shifman, Z. Phys. C9 (1981) 347; P. Pascual and R. Tarrach, 
Phys. Lett. B113 (1982) 495. 

J. Sexton, A. Vaccarino and D. Weingarten, Phys. Rev. Lett. 75 (1995) 
4563. 

Y. Kim, H.K. Lee and M. Rho, Phys. Rev. C52 (1995) R1184. 



S.M. Johns, P.J. Ellis and J.M. Lattimer, |nucl-th/ 960400^ , Astrophys. 
J. 473 (1996) in press. 



[15] M. Abramowitz and LA. Stegun, eds. Handbook of Mathematical Func- 
tions (Dover, NY, 1965). 



[16] T.D. Lee and M. Margulies, Phys. Rev. Dll (1975) 1591. 

[17] P. Gerber and H. Leutwyler, Nucl. Phys. B321 (1989) 387. 

[18] F.R. Brown, P.P. Butler, H. Chen, N.H. Christ, Z. Dong, W. Schaffer, 
L.I. linger and A. Vaccarino, Phys. Rev . Lett. 65 (1990) 2491. 



[19] H. Satz, hep-ph/ 96 113661 , Proceedings of the First Catania Relativistic 



Ion Studies (Acicastello, Italy, 1996), to be published. 
[20] N.O. Agasyan, JETP Lett. 57 (1993) 208. 

Figure Captions 

Figure 1. The mean sigma field, v = a/ao, as a function of temperature. 
The dashed (solid) line corresponds to the presence (absence) of explicit 
chiral symmetry breaking. The dotted curve indicates a thermodynamically 
unstable region. 

Figure 2. The mean glueball field, x = 0/00, as a function of temperature. 
The dotted and dash-dotted curves indicate thermodynamically unstable re- 
gions. See caption to Fig. 1. 

Figure 3. The pion effective mass as a function of temperature. See caption 
to Fig. 1. 

Figure 4. The sigma effective mass as a function of temperature. See 
caption to Fig. 1. 

Figure 5. Schematic representation of the grand potential per unit volume, 
^/V, as a function of the glueball field x fo^' various temperatures with 
e[ > 0. Note that for each case Q is normalized to zero at x = 0- 
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Figure 6. The glueball effective mass as a function of temperature. See 
caption to Fig. 1. 

Figure 7. Plot of the pressure versus energy density for the temperatures 
indicated. See caption to Fig. 1. 

Figure 8. Representative diagram whose thermal part is evaluated to de- 
termine {{tp'^y). 
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